BMC Biology

Conserved synteny enables the reconstruction of ancestral chromosome organization across deeply diverged lineages, as demonstrated across metazoans. In Spiralia, comparative genomic analyses have suggested an ancestral karyotype comprising 20 chromosomes. We report high-quality genome assemblies for two Lima species (Lima caribaea and Lima lima) with 20 chromosomes and show that their chromosome organization matches this ancestral state, including four previously described fusion-with-mixing events. Interestingly, other Limidae genus did not retain the ancestral karyotype: they have undergone subsequent lineage-specific fusions. The conserved organization of Lima chromosomes provides a rare genomic anchor for reconstructing ancestral chromosome organization and lineage-specific chromosome evolution across Spiralia and, more broadly, Metazoa.

The twisted-wing parasites (order Strepsiptera) are well-known for their unusual life histories as obligate endoparasites of a wide variety of other insect orders. Here, we present the first high-quality genome assembly of a leafhopper-parasitizing strepsipteran species. We sequenced the genome of Elenchus koebelei using PacBio HiFi technology, yielding a primary assembly of 63.7Mb in 96 contigs, with an N50 of 1.60Mb and a BUSCO completeness score of 82.3%. In total, we recovered 8,167 protein-coding gene annotations, and 34.2% of the genome annotated as repeats. For the first time, we compare the genome content of three strepsipteran species and find that they are remarkably consistent in broad structure, with approximately 600 insect universal single-copy orthologs that appear to be completely missing from the order. These similarities imply a conserved genomic blueprint for strepsipteran parasitic strategy and/or differences that must be investigated at the post-translational level.

Amoebic gill disease (AGD), caused by the parasite Neoparamoeba perurans, is a major welfare and economic burden for the Atlantic salmon aquaculture industry. However, the spatial organisation of immune functions within the gill, the site of infection, remains poorly resolved. Here, we present the first spatial transcriptomic atlas of the Atlantic salmon gill affected by AGD, generated using 10x Genomics Visium on four AGD-positive fish, confirmed by gill scoring and quantitative PCR for parasite load. As no uninfected gills were included, this study maps the anatomical and immunological organisation of AGD-affected gills rather than defining AGD-specific pathology. Unsupervised clustering of spatial gene expression resolved distinct transcriptional domains that grouped into four principal anatomical compartments of the gill, including (i) epithelial and lamellar regions defined by keratin and osmoregulatory gene expression, (ii) gill arch regions enriched for connective tissue and smooth muscle markers, (iii) an erythroid-associated vascular cluster, and (iv) a skeletal muscle cluster with signatures of glycolysis and oxidative phosphorylation. Across these domains, we identified a spatially organised and diffused immune landscape, with pattern recognition receptors, antimicrobial effectors and chemokines concentrated at the filament surface, and immunoglobulin light chain transcripts consistent with B lymphocyte activity enriched toward the gill arch. Integration of the spatial data with a published single-nucleus reference atlas using reference-based deconvolution resolved the cellular identity underlying these domains and revealed frequent co-localisation of immune and structural cell types, suggesting regulated cell to cell communication within the gill-associated lymphoid tissue. Overall, our results show that the Atlantic salmon gill is not a transcriptionally uniform surface but a regionally specialised organ, with immune functions distributed across, and specialised within, distinct anatomical compartments

The origin of multicellular animals from unicellular ancestors remains a major evolutionary question. To address this transition, a reconstruction of the nature of the last unicellular ancestor of animals (LUAA) is key, and requires broad sampling of extant unicellular relatives, which together with animals form the clade Holozoa. Environmental surveys suggest that unicellular Holozoa diversity extends beyond the five groups with characterized members (Choanoflagellatea, Ichthyosporea, Filasterea, Corallochytrea, and Atreyea). However, the extent of the environmental diversity within and outside these groups remains unclear. To fill this gap, we here systematically re-evaluated Holozoa diversity using environmental 18S rDNA sequences compiled from updated datasets, including long-read metabarcoding surveys. We report nine uncharacterized lineages inside Holozoa (seven newly identified here), provisionally placed outside the described groups. We also uncovered a previously unreported diverse lineage deeper inside Opisthokonta, separate from Holozoa and Fungi. In addition, we identified 21 uncharacterized lineages within Choanoflagellatea and Ichthyosporea (eleven newly reported here). We used short-read metabarcoding datasets to characterize the environmental and geographical distributions of these lineages. Many are globally distributed and some can reach relative abundances comparable to those of the described groups, while others are restricted to specific environments. Our findings reveal substantial uncharacterized diversity missing from the cultured taxa used in comparative studies of animal origins. By placing this diversity in a phylogenetic and environmental context, our study will guide the targeted isolation efforts needed to broaden the taxon sampling available to reconstruct LUAA and unravel the unicellular-to-multicellular transition that gave rise to animals.

RNA interference (RNAi) offers a sequence-specific approach to pest control. In insects, Dicer-2 processes double-stranded RNA (dsRNA) into small interfering RNA (siRNA) duplexes, from which the RNA-induced silencing complex (RISC) retains a guide strand. Only antisense-loaded RISC can mediate cleavage of the target transcript. However, how sequence features shape the RISC-bound siRNA pool in pests remains poorly understood, limiting opportunities for sequence optimization. Here, we profiled RISC-bound siRNAs following injection of 34 insecticidal dsRNAs targeting 11 essential genes in Tribolium castaneum larvae. We computationally reconstructed 7,879 siRNA pairs and examined associations between sequence features and strand bias. Differences in GC identity at terminal paired positions 1 to 5, used as a proxy for local thermodynamic asymmetry, correlated with strand bias, with the strongest correlations at the first two paired positions. ORF targeting and reduced predicted antisense self-folding were also associated with higher antisense fractions. Analysis of non-templated terminal additions revealed predominantly 3-prime uridylation, a known signature of small RNA turnover, along with putative 3-prime trimming. Among ORF-associated siRNA pairs, sense strands showed higher relative U-tailing abundance, based on 3-prime uridylated and putatively trimmed-and-3-prime-uridylated reads relative to perfect 21-nt reads, than antisense strands. Antisense strands with the least predicted self-folding also showed low relative U-tailing abundance. These observations are consistent with sequence-dependent contributions from both guide-strand selection and differential post-RISC-loading siRNA retention, although a causal link remains to be established. The identified associations provide a basis for testing whether dsRNA sequence optimization can improve pest control efficacy and reduce off-target activity.

Adult stem cells can survive diverse insults, from ionizing radiation to toxic chemicals. In vertebrates, ABC transporter expression and dye-efflux capacity have been associated with stem-cell protection and multidrug resistance. How compound accumulation differs between adult stem cells and their differentiated daughters within an intact tissue remains poorly understood. Here, we show that adult Drosophila intestinal stem/progenitor cells maintain a distinct low-accumulation state for membrane-permeable compounds and that ABC transporters contribute to this phenotype and to the response to cytotoxic stress. Across nine membrane-permeable probes, stem/progenitor cells accumulated eight at lower levels than differentiated enterocytes. We selected tetramethylrhodamine methyl ester (TMRM), a membrane-permeable fluorescent probe that showed a strong difference in accumulation between stem/progenitor cells and enterocytes, as a quantitative readout for an RNAi screen of all 56 Drosophila ABC-family genes. The screen identified seven candidate transporters whose depletion increased the stem/progenitor-to-enterocyte TMRM fluorescence ratio. Further characterization of one of them, ABCB7, showed that ABCB7-depleted guts contained fewer mitotic cells during cytotoxic drug exposure. These results implicate ABC transporters in lower compound accumulation in Drosophila intestinal stem/progenitor cells and support a role for ABC transporters in their response to cytotoxic challenge.

Background Comparisons of single-cell and single-nucleus RNA sequencing (scRNA-seq and snRNA-seq) data have been described in some mammalian tissues and, subsequently, in Drosophila, but remain unexplored in most invertebrate lineages. The xenacoelomorphs occupy key phylogenetic positions, yet they differ anatomically from mammals. They have a reduced extracellular matrix, high-salt body fluid, and no circulatory system. Despite these differences, they possess a well-developed nervous system. One of the xenacoelomorphs, the photosymbiotic acoel (Praesagittifera naikaiensis) also harbours symbiotic Tetraselmis algae, whose RNA can be co-captured with host RNA. Results We compared scRNA-seq and snRNA-seq data from whole P. naikaiensis specimens. Both methods yielded high-quality data with comparable gene detection but a larger share of scRNA-seq reads derived from symbiotic algae. Gene-level analyses revealed that neural genes were enriched in snRNA-seq relative to non-neural genes. Cross-method label transfer and integration-based validation identified six snRNA-seq clusters, including some neural populations, that lacked a clear scRNA-seq counterpart. In contrast, three cell populations, including muscle and metabolically active clusters, were reciprocally validated as captured by both methods. Conclusions Our results show that key snRNA-seq advantages, particularly the enhanced recovery of neural transcripts, are recapitulated in our dataset, which is consistent with previous reports in mammals. Furthermore, snRNA-seq reduces symbiont-derived reads and recovers several cell populations underrepresented in scRNA-seq. These findings provide practical guidance for cell atlas construction in non-model, symbiotic invertebrates.

Prospero-related homeobox (Prox) genes encode a family of transcription factors that play essential roles in the development of several organs and systems, including the central nervous system, lymphatic endothelium, musculature, and liver. Despite their developmental importance, the evolutionary history of the vertebrate Prox gene family remains poorly understood. In this study we combined phylogenetic and synteny analysis to characterise the evolution of the Prox family in vertebrates. Our results reveal that two to three Prox subfamilies were already present in the last common ancestor of jawed vertebrates. We clarify the identity and evolutionary relationships of well-studied members of this family and identify multiple independent losses of Prox2 and Prox3 genes in specific vertebrate lineages. Furthermore, we uncover evidence for the existence of a fourth Prox gene in the ancestral vertebrate genome, which was subsequently lost. Overall, this study provides the first comprehensive analysis of the evolutionary history of the vertebrate Prox gene family and establishes a foundations for future studies on the functional roles of these genes.

Double-stranded RNA (dsRNA) biopesticides that silence essential fungal genes by spray-induced gene silencing (SIGS) are advancing towards registration, underpinned by the claim that silencing is confined to sequence-matched target organisms, a claim that has never been tested systematically against beneficial fungi. We screened twelve dsRNA constructs (the whole transcripts of eleven canonical antifungal target genes of Fusarium graminearum: three CYP51 sterol 14-demethylases, six chitin synthases and two {beta}-tubulins, plus a concatenation of the three CYP51 transcripts modelling the flagship CYP3RNA design) for perfect 21-nucleotide (21-nt) identity, the primary match criterion of recent regulatory bioinformatics frameworks, against the transcriptomes of seven non-target fungi, including commercial biocontrol agents (Trichoderma harzianum, T. virens, Beauveria bassiana, Metarhizium anisopliae, M. brunneum), the arbuscular mycorrhizal fungus Rhizophagus irregularis and Saccharomyces cerevisiae. Eleven of the twelve constructs carried off-target hazard; only CYP51C was completely clean. TUB2b carried 580 perfect 21-nt matches and contained no hazard-free window of 250 consecutive 21-mer start positions, whereas 26.6-60.0% of screened windows were hazard-free in the designable genes. Hits landed overwhelmingly on orthologues of the targeted gene and concentrated in the Hypocreales; R. irregularis and S. cerevisiae were nearly clean. These are hazard findings, not risk findings: perfect siRNA-length matches to beneficial fungi are common in canonical antifungal dsRNA targets, and current regulatory screening panels contain no fungus that would detect them.

Background: Nervous system regeneration requires creation of new neurons and their integration into new and existing networks. How new neurons find their targets and reconnect in the environment of regenerating tissue remains a major gap in knowledge. The planarian Schmidtea mediterranea is well suited to the study of neural regeneration due to the remarkable ability of planarians to regenerate their entire body, including the nervous system. However, we have not yet been able to explore synaptogenesis in regeneration using planarians due to a lack of tools to study synapses. Here we investigate potential players in planarian synapse, focusing on the conserved syntaxin protein family, which functions in neurotransmitter vesicle fusion at the synapse. Results: We found that reduction of either syntaxin-1 (stx1) or syntaxin-binding protein 1-1 (stxbp1-1) decreased locomotion and impeded brain regeneration, specifically in synapse-rich regions. Conclusions: Our observations support STX-1 and STXBP1-1 as bona fide synapse components, providing a starting point for studies neuronal connectivity during homeostatic remodeling and regeneration.

Early embryos must mount stress defenses without compromising the developmental programs that drive morphogenesis, but whether genetic background alters this balance remains unclear. We compared cellularizing Drosophila melanogaster embryos from inbred and outcrossed backgrounds under control conditions and after acute or chronic heat-stress. Both backgrounds developed successfully and showed modest transcriptomic differences at 18{degrees}C. Heat stress at 32{degrees}C, however, sharply separated their phenotypic and transcriptional responses. Inbred embryos showed poor larval hatching after acute and chronic heat, and this deficit persisted through adult eclosion. Their transcriptomes exhibited a greater stress-response burden, with enrichment of proteostasis, DNA repair, mitochondrial maintenance, RNA processing, and metabolic remodeling programs. Outcrossed embryos retained substantially greater developmental survival and higher abundance of transcripts associated with development, adhesion, epithelial organization, signaling, transport, sensory responses, and morphogenesis. Background-dependent divergence was greatest during acute heat. Under chronic heat, background-dependent differences were fewer, with most overlapping those seen during acute heat and a smaller set detected only under chronic heat. Transcript-origin analysis further showed that acute outcrossed-biased genes were predominantly maternally derived and showed no enrichment for purely zygotic transcripts, whereas chronic outcrossed-biased genes contained a larger zygotic component. In summary, the less viable inbred background mounted the broader molecular stress signature, whereas the more viable outcrossed background retained greater abundance of development-associated transcripts. These findings identify genetic background as a major determinant of how embryonic heat responses are distributed between molecular stress defense and developmental gene expression.

Two rounds of whole-genome duplication (2R-WGD) in the vertebrate lineage are thought to have driven gene diversification. Mammals have four genes encoding five nebulin superfamily proteins (LASP1, LASP2/NEBL, NEB, and NRAP) hypothesised to have arisen from invertebrate LASP via WGD. Here, we examined the evolutionary fate of the nebulin superfamily after an additional WGD, teleost-specific genome duplication (TGD). Genome and transcriptome analysis of zebrafish (Danio rerio) and medaka (Oryzias latipes), of which the complete genomes were recently determined, showed only lasp family genes retained duplicate copies. Zebrafish retained two LASP2/NEBL-derived ohnologs; one (LASP2L) had been misidentified as a nebulette gene in the Zebrafish Information Network (ZFIN) despite lacking nebulette-specific exons. Another (LASP2/NEBL) is newly identified in this study because it lies in a newly determined region of the latest genome assembly. Contrastingly, medaka retained two LASP1 ohnologs and one LASP2/NEBL. The duplicated lasp family proteins showed overlapping and specific tissue expression patterns. We confirmed all three zebrafish lasp proteins bound actin and co-localized with actin-rich structures when ectopically expressed in cultured fibroblasts. This supports that TGD-derived lasp duplicates underwent lineage-specific retention and expression divergence in teleosts whilst conserving core actin-binding properties.

The olfactory ventricle (OV), an anterior extension of the lateral ventricular system into the olfactory bulbs, is generally considered a transient developmental structure that closes during early postnatal life. However, this view largely derives from a limited number of laboratory models, and the evolutionary distribution and cellular organization of the adult OV remain poorly understood. Here, we combined a comparative survey of mammalian neuroanatomy with phylogenetic analyses and cellular characterization of the OV across species. We compiled anatomical data for 154 mammalian species spanning 26 orders and 80 families, revealing widespread persistence of a patent OV into adulthood. Phylogenetic analyses supported OV presence as the ancestral mammalian state and indicated that its absence is phylogenetically structured and consistent with repeated secondary losses. Chiroptera displayed multiple lineage-specific losses, whereas variation within Muridae extended to the intraspecific level, with contrasting OV persistence among Mus musculus and Rattus norvegicus populations. Comparative histological analyses in sheep, Peromyscus maniculatus, and CD1 mice revealed that, despite marked interspecific differences in OV morphology, the persistent ventricular region retains a conserved cellular organization characterized by progenitor- and glial-associated markers, proliferative cells, immature neuronal populations, and a differentiated apical/ciliary interface. Although these observations do not establish local neuronal production or demonstrate that the OV constitutes an autonomous neurogenic niche, their convergence across species indicates that the persistent adult OV is not simply an anatomically residual lumen but a specialized component of mammalian olfactory bulb organization. Together, our findings identify the OV as an ancestral and widely conserved feature of the mammalian brain and provide an evolutionary and cellular framework for investigating the developmental mechanisms controlling its persistence and its potential contribution to adult olfactory neurogenesis and plasticity.

Complex traits that span multiple tissues and systems often integrate large numbers of genes across development, physiology, and behavior, making it challenging to identify their essential components. Blood feeding in mosquitoes is one such trait. It is ancestral to the mosquito family, maintained in most species for ~200 million years, and was independently lost in three lineages. These convergent losses offer a natural experiment to discover the genetic, physiological, and neural features required for blood feeding. We assembled high-quality, chromosome-level genomes for seven mosquito species, along with whole-brain tomographic reconstructions. Our study spanned the three known non-blood-feeding lineages (Toxorhynchites rutilus, Topomyia yanbarensis, and Malaya genurostris), blood-feeding relatives, and the variable blood feeder Wyeomyia smithii. Comparing orthologous gene clades, we detected convergent gene loss specific to the three lineages that had lost blood feeding. The losses include the salivary platelet-aggregation inhibitor Aegyptin, blood-activated serine proteases such as Chymotrypsin-1 and 2, and a carboxylesterase expressed in the female fat body and brain glia. The loss of blood feeding also extended to chemosensation. Non-blood feeders lack two odorant-binding protein clades, two ionotropic receptor clades associated with blood-component taste detection, and odorant receptor clades expressed in a discrete, strongly female-biased population of antennal neurons. Female-biased head gene expression was reduced in non-blood feeders. Finally, examination of whole-brain tomographic reconstructions across the species revealed smaller antennal lobes in non-blood-feeding females, consistent with reduced olfactory input. Together, these findings identify a compact set of genes, expression patterns, and brain regions associated with blood feeding, offering an evolutionary entry point for functional dissection of how this complex and dangerous trait is built and dismantled.

Lateral gene transfer (LGT) from organellar to eukaryotic genomes is ubiquitous, resulting in the presence of widespread nuclear mitochondrial DNA segments (NUMTs) and/or plastid gene transfers. While LGT from bacterial associates of eukaryotes is generally less common, many arthropod genomes are littered with LGT originating from Wolbachia (order Rickettsiales) - a vertically transmitted, obligate intracellular symbiont of invertebrates. This contrasts with Rickettsia, a related genus including many important human pathogens, which has not been shown to promulgate LGT within its arthropod vectors. Here, we present the 8.6 Gb genome of a cell line derived from A. variegatum (the tropical bont tick), which is the primary vector of Rickettsia africae (agent of African tick-bite fever). In addition to numerous NUMTs and diverse families of transposable elements, an almost-complete chromosome of R. africae origin was identified in the cell line genome, localised to the putative sex chromosome. Sequencing of field-collected A. variegatum confirmed the presence of this large LGT in wild ticks, although the absence of transferred genes from the R. africae plasmid provided a means to differentiate LGT from genuine rickettsial infections. These findings highlight the imperative to consider the possibility of LGT when screening vectors for pathogens by PCR.

Clownfishes are a complex of 28 species that rapidly diversified after acquiring mutualism with host sea anemones around 15 million years ago. The genomic mechanisms behind this radiation are only beginning to emerge, and structural variants (SVs) have received little attention. Previous phylogenomic studies suggested two large inversions on chromosome 18, but their distribution across the clownfish phylogeny and their potential effects have not been formally explored. SVs can disrupt gene function and regulation or alter gene dosage, and they have been shown to play a central role in adaptive evolution and diversification across many taxa. Here, we characterized SVs across the clownfish radiation using long-read PacBio sequencing, generating high-quality genome assemblies for sixteen species spanning the main lineages of the genus. We identified and catalogued nearly 130,000 SVs spanning close to 570 Mb. Short insertions and deletions were the most numerous, but large inversions accounted for most of the affected sequence. No SV was consistently associated with host specialization. We nevertheless confirmed a single inversion of approximately 17 Mb on chromosome 18, clarifying earlier reports of two separate rearrangements at this locus, and identified an additional inversion on chromosome 9 whose distribution suggests a history of hybridization. This study provides the first genome-wide characterization of structural variation in clownfishes, revealing extensive SVs despite the group's rapid and recent radiation. It also delivers new, high-quality genomic resources that will support future research into this iconic group of coral reef fishes and their adaptation to sea anemones.

Host specialization in herbivorous insects is often associated with divergence in chemosensory abilities. Here, we investigated the possible contribution of odorant receptors (ORs) in host plant restriction in the lily moth Spodoptera picta, a species specialized on Amaryllidaceae. Manual annotation of S. picta ORs in its genome revealed a repertoire similar in size and composition to those of its polyphagous sister species, S. littoralis and S. litura, suggesting that specialization did not involve major gene loss or expansion in the lily moth. To assess functional divergence beyond gene number, we applied a large scaled structure-based virtual screening approach to the entire OR repertoires of these three Spodoptera species, generating ligand-binding profiles for 120,591 volatile compounds. Among 69 1:1:1 OR orthologs, 24 exhibited divergent predicted binding spectra. We pinpointed OR29 that we also found to be highly expressed in both male and female antennae of S. picta through a RNAseq approach. Functional assays demonstrated that S. picta OR29 acquired heightened sensitivity to limonene enantiomers, volatiles emitted by host Amaryllidaceae inflorescences. Site-directed mutagenesis revealed that a single amino acid substitution within the predicted binding region underlies this shift in sensitivity. These results show that host specialization in S. picta has not been accompanied by significant OR repertoire remodeling, but rather by subtle molecular changes that fine-tune receptor sensitivity to host-derived volatiles.

Menopause timing links reproductive lifespan with later life health, motivating investigation of its molecular and hematological correlates. We combined Generation Scotland analyses of women sampled at or after reported menopause (n = 1,444) and initially menstruating women followed for GP ascertained menopause onset (n = 1,725; 251 events) with bidirectional and multivariable Mendelian randomization. Higher age adjusted GrimAge V1, GrimAge V2 and DunedinPACE corresponded to 0.34 to 0.40 years earlier menopause per SD in the retrospective cohort at nominal P < 0.05. Prospective associations were nonsignificant, and cohort and primary bidirectional aging-related genetic associations did not survive multiple-testing correction. Genetically predicted later menopause was associated with higher neutrophil percentage and count and lower lymphocyte percentage in women; these associations and female vs. male contrasts survived correction. CRP conditional on menopause timing was associated with higher GrimAge and PhenoAge acceleration and shorter leukocyte telomere length after correction. None of 36 estimated inflammatory indirect effects reached P < 0.05. These findings highlight blood-cell composition as a prominent correlate of reproductive timing and identify selected methylation-aging signatures for further study. Together, they position reproductive hematological relationships as a specific avenue for understanding the connection between menopause and systemic aging.

Torpor is a physiological strategy used by which some species conserve energy by lowering metabolic rates and body temperatures. While mammalian hibernation is well-studied, the molecular mechanisms governing the rapid entry into avian torpor remain poorly understood. This study represents the first multi-tissue whole-transcriptomic analysis of torpor entry in an avian system. By combining real-time thermal imaging and respirometry with transcriptomic profiling across seven tissues, we identified the key molecular pathways orchestrating this transition in Anna's hummingbirds (Calypte anna). Our results reveal that avian torpor entry is characterised by a subtle transcriptomic change (<5% of the genome differentially expressed across tissues) compared to mammalian hibernation (10%-60%). Most of these genes are differentially expressed before torpor entry is complete, rather than tracking body temperature linearly. This indicates that active transcriptional reprogramming precedes and potentially drives the transition. Thus, hummingbird torpor entry is not a passive shutdown, but a highly regulated, tissue-specific process. Vital organs, such as the heart and lungs, remain transcriptionally stable to maintain essential functions, while the liver and gut exhibit extensive metabolic rewiring. We found evidence of global transcriptional and translational suppression, and cell-cycle arrest: the three major contributors to cellular energy budgets. We hypothesise that avian torpor entry is orchestrated through an arrest of apoptosis, a metabolic switch from carbohydrate to lipid utilisation, coupled with alterations in alternative splicing and circadian rhythm regulation. Finally, the counterintuitive upregulation of genes involved in mitochondrial metabolism suggests a preparation for rapid arousal despite the depressed metabolic state.

Whole-genome duplications (WGDs) have long been proposed as evolutionary facilitators of biological complexity and diversification. Ancient WGDs occurred deep in the evolutionary history of both plants and vertebrates, but appear to be rare in invertebrates. However, it is unclear whether this reflects genuine scarcity or inadequate detection methodologies, and current methods may be insufficient to detect very old WGDs where extensive gene loss has followed. Here, we used a combination of paralogous and orthologous synteny-based methods incorporating bilaterian ancestral linkage groups (ALGs) to search for genomic signatures of WGD in members of the invertebrate phylum Bryozoa. This revealed strong evidence for an ancient WGD in the freshwater bryozoan Cristatella mucedo, Class Phylactolaemata, with duplicated paralogous regions across every chromosome. Only ~10% of duplicated genes are retained, including a Hox cluster duplication. All eight genomes from the bryozoan class Gymnolaemata also show signatures of ancient WGD in the form of 2:1 ratios in orthologous synteny comparisons, which remain detectable despite only ~5% duplicate retention. Gene trees support a scenario in which all extant bryozoans share an ancient WGD and much of the genome rediploidised independently in the two clades. Duplicates retained after bryozoan WGD are enriched in cilia genes and are preferentially expressed in the lophophore, the tentacular feeding organ covered in highly specialised cilia. We argue that cryptic WGD in bryozoan evolution facilitated adaptation to a sessile, filter-feeding lifestyle, and propose that incorporation of ALGs adds to the power of synteny-based WGD-detection methods.